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Microwave spectroscopy was used to probe the superfluid-Mott Insulator tran- 
sition of a Bose-Einstein condensate in a 3D optical lattice. Using density de- 
pendent transition frequency shifts we were able to spectroscopically distin- 
guish sites with different occupation numbers, and to directly image sites with 
occupation number n = 1 to n = 5 revealing the shell structure of the Mott 
Insulator phase. We use this spectroscopy to determine the onsite interaction 
and lifetime for individual shells. 

The Mott insulator transition is a paradigm of condensed matter physics. It describes how 
electron correlations can lead to insulating behavior even for partially filled conduction bands. 
However, this behavior requires a commensurable ratio between electrons and sites. If this con- 
dition on the density is not exactly fulfilled, the system will still be conductive. For neutral 
bosonic particles, the equivalent phenomenon is the transition from a superfluid to an insula- 
tor for commensurable densities. In inhomogeneous systems, as in atom traps, the condition 
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of commensurability no longer applies: for sufficiently strong interparticle interactions, it is 
predicted that the system should separate into Mott insulator shells with different occupation 
number, separated by thin superfluid layers EIEIEI. 

The recent observation of the superfluid to Mott Insulator (MI) transition with ultracold 
atoms has stimulated a large number of theoretical and experimental studies [see Q and 
references therein]. Atomic systems allow for a full range of control of the experimental param- 
eters including tunability of the interactions and defect-free preparation making them attractive 
systems for studying condensed matter phenomena. The Mott insulator phase in ultracold atoms 
has been characterized by studies of coherence, the excitation spectrum, and noise correlations 
||4l|6l|7|. Recently, using spin-changing collisions, lattice sites with two atoms were selectively 
addressed and the suppression of number fluctuations was observed |l8l|. 

The most dramatic new feature of the MI phase in ultracold atoms is the layered structure 
of the Mott shells, however until now this feature has not been directly observed. In this paper, 
we combine atoms in the Mott insulator phase with the high resolution spectroscopy used for 
atomic clocks, and use density dependent transition frequency shifts to spectroscopically resolve 
shells with occupancies from n = 1 to 5 and directly image their spatial distributions. 

Bosons with repulsive interactions in an optical lattice are accurately described by the 
Hamiltonian ED, 

H = - JY^ a\aj + ]:UYni{hi - 1) + ^(ei - /i)^,, (1) 

where the first two terms are the usual Hamiltonian for the Bose-Hubbard model and the last 
term adds in the external trapping potential, and where J is the tunneling term between nearest 
neighbors, a\ and are the boson creation and destruction operators at a given lattice site, 
w{x)\'^d^x is the repulsive onsite interaction, where m is the atomic mass, 
a is the s-wave scattering length and w{x) is the single particle Wannier function localized to 
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the i*" lattice site, and = a-di is the number operator for bosons at site i. The last term in the 
hamiltonian is due to the external trapping confinement of the atoms, where = Vext{ri) is the 
energy offset at the ith site due to the external confinement and ji is the chemical potential. 

The behavior of this system is determined by the ratio of J/U. For low lattice depths, the 
ratio is large and the system is superfluid. For larger lattice depths, the repulsive onsite energy 
begins to dominate, and the system undergoes a quantum phase transition to a MI phase. For 
deep lattices the atoms are localized to individual lattice sites with integer filling factor n. This 
filling factor varies locally depending on the local chemical potential /Xj = — as 

n = Modii^i/U), (2) 

and decreases from the center to the edge of the trap. 

To prepare the atoms in the Mott insulating phase, a *^Rb Bose-Einstein condensate was first 
created in the the \F — 1, = —1) state using a combination of a loffe-Pritchard magnetic 
trap and an optical dipole trap. The optical trap was oriented perpendicular to the long axis 
of the magnetic trap, creating a more isotropic trapping potential which was better matched 
to the optical lattice. The laser beam for the optical trap had a 1/e^ waist of 70 iim, and 
was retro-reflected. However, the polarization of the retro-reflected beam was rotated such that 
the interference between the two beams had minimal contrast. The resulting trap had radial 
(axial) trap frequencies of u; = Inx 70 (20) Hz, where the axial direction is now parallel to 
the optical trap. A 3D optical lattice was created by adding two additional retro-reflected laser 
beams derived from the same laser at A = 1064 nm. The lattice was adiabatically ramped up by 
rotating the polarization of the retro-reflected optical trapping beam to increase the interference 
contrast along that axis and by increasing the laser power in the other two axes. The lattice depth 
was increased using an exponential ramp with a 40 ms time constant. All three beams were 
linearly polarized orthogonal to each other and had different frequency detunings generated 
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using acousto-optic modulators. The lattice depth was up to 40 Erec, where Erec = h'^k'^/ (2m) 
and /c = 27r/A is the wavevector of the lattice. At 40 Erec the lattice trap frequency at each site 
was cuiat = 27r X 25 kHz, and the external trap frequencies increased to = 27r x 110 (30) Hz in 
the radial (axial) direction. 

Zeeman shifts and broadening of the clock transition from the F = 1 to the F = 2 state 
were avoided by using a two photon transition between the |1, — 1) and the |2, 1) state, where at 
a magnetic bias field of ~ 3.23 G both states have the same first order Zeeman shift lITOll . The 
two photon pulse was composed of one microwave photon at a fixed frequency of 6.83 GHz, 
and one rf photon at a frequency of around 1 .67 MHz. The pulse had a duration of 100 ms, and 
when on resonance the fraction of atoms transferred to the |2, 1) state was less than 20%. After 
the pulse, atoms in the |2, 1) state were selectively detected with absorption imaging using light 
resonant with the 5^51/212, 1) 5^P3/2|3, 1) transition. For observing the spatial distribution 
of the Mott shells, the atoms were imaged in trap. For recording spectra, the atoms were released 
from the trap and imaged after 3 ms of ballistic expansion in order to reduce the column density 

When the two photon spectroscopy is performed on a trapped condensate without a lattice, 
the atoms transferred to the |2, 1) state have a slightly different mean field energy due to the 
difference between the a2i and an scattering lengths, where 021 is the scattering length between 
two atoms in states |2, 1) and |1, —1), and an is the scattering length between two atoms in the 
|1, — 1) state. This difference in the scattering lengths leads to a density dependent shift to the 
resonance frequency oc p(a2i —o-u), where p is the condensate density ifTOl . This coUisional 
shift is commonly referred to as the 'clock' shift ifTTl due to its importance in atomic clocks, 
where cold collisions currently limit the accuracy llT2llT3l . When performed on a condensate 
with peak density po in a harmonic trap, in the limit of weak excitation, the line shape for the 
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2-photon resonance is given by HH : 



Ap,AE f-^^' 

with the mean field energy difference 

AE= (a2i - an) • (4) 

rem 

In the case of ^^Rb, 021(011) = 5.19 (5.32) nm ifTSl . Both the frequency shift and the linewidth 
increase with the condensate density. As the lattice is ramped on, the peak density of the 
condensate in a given lattice site increases as 

Po{r) = (^/i - ^mul^pr^^ 1/U, (5) 

where ujtrap is the external trap frequency for the combined magnetic and optical trap, and using 
the Thomas Fermi approximation /i, the chemical potential is given by 

"^U vii j ■ 

where is the total atom number. For low lattice depths the system is still a superfluid, delo- 
calized over the entire lattice. However, the two-photon resonance line is shifted and broadened 
due to the increased density, with the center of the resonance at v = + 2poAE /3h. For deep 
lattices in the MI regime the repulsive onsite interaction dominates, number fluctuations are 
suppressed, and each lattice site has a sharp resonance frequency determined by the occupation 
number in the site. The separation between the resonance frequencies for the n and n — 1 MI 
phases is given by 

= ^{a2i- ail) /au. (V) 

The linewidth of the resonances is no longer broadened by the inhomogeneous density and 
should be limited only by the bandwidth of the two photon pulse. 



Fig. [T] shows the transition from a broadened resonance to several sharp lines as the lattice 
depth was increased. At a lattice depth of V = 5 Erec, the line was broadened and the line 
center was shifted slightly due to the increased density. At V = 10 Erec, the line was shifted 
and broadened further, and in addition the line shape became asymmetric as the atom number 
in lattice sites with small occupation was squeezed. For deeper lattice depths the system under- 
went a phase transition to a MI phase and discrete peaks appeared corresponding to MI phases 
with different filling factors, and for V = 35Erec, MI phases with occupancies of up to 5 were 
observed. 

When the lattice depth was increased inside the MI regime (from V = 25 Erec to 35 Erec), 
the separation between the resonance peaks increased presumably due to the larger onsite inter- 
action energy as the lattice trap was increased. As given in Eq. (7), the separation between the 
peaks provides a direct measurement of the onsite interaction energy U. As shown in Fig. |2^, 
our results are in good agreement with calculated values of U. Fig. |2b shows that although 
the separation between the n = 1 , 2 and 3 peaks is approximately constant, for higher filling 
factors the separation between the peaks decreases; the effective onsite interaction energy be- 
comes smaller for higher filling factors. This shows that for low occupation numbers the atoms 
occupy the ground state wavefunction of the lattice site, whereas for larger occupation numbers, 
the repulsive onsite interaction causes the wavefunction to spread out lowering the interaction 
energy. From a variational calculation of the wavefunction similar to llT6l . we find the onsite 
energy for the n = 5 shell should be ~ 20% smaller than that for the n = 1 shell, is in agreement 
with the measured value shown in Fig.|2b. 

The peaks for the different occupation numbers were spectrally well separated. Therefore, 
on resonance, only atoms from a single shell were transferred to the |2, 1) state. An image of 
these atoms (without any time-of-flight) shows the spatial distribution of this shell. Fig. [Sj), 
shows absorption images for n = 1 - 5 shells. As predicted lU, the n = 1 MI phase appears 
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near the outer edge of the cloud. For larger n, the radius of the shell decreases, and the n = 5 
sites form a core in the center of the cloud. The expected radius for each shell was obtained 
from Eq. 2 using the measured values for the onsite interaction. The observed radii were in 
good agreement, except for the n = 1 shell which may have been affected by anharmonicities in 
the external trap. Absorption images taken with rf frequencies between the peaks show a small 
signal, which may reflect the predicted thin superfluid layers between the insulating shells, 
however this needs be studied further. The expected absorption image of a shell should show a 
column density with a flat distribution in the center and raised edges. However due to limitations 
(resolution and residual fringes) in our imaging system, these edges were not resolved. 

Since we were able to address the different MI phases separately, the lifetime for each shell 
could be determined. For this, the atoms were first held in the lattice for a variable time r before 
applying the 100 ms two-photon pulse. For the n = I MI phase, ignoring technical noise, the 
lifetime should only be limited by spontaneous scattering from the lattice beams. Even for the 
deepest lattices, the spontaneous scattering rate is < 10^^ Hz. For the n = 2 MI phase the 
lifetime is limited by dipolar relaxation, which for ^^Rb is slow with a rate < 10^^ Hz. For 
sites with n > 3 the lifetime is limited by 3-body recombination with a rate 'jn(n — l){n — 2) 
ifTTll with 7 = 0.026 Hz for our parameters. This gives 3-body lifetimes of of t^b = 6.2 s, 1.6 s, 
and 0.6 s for the n = 3, 4, and 5 MI phases respectively. This calculation of 7 assumes for the 
density distribution the ground state of the harmonic oscillator potential, so for higher filling 
factors the actual lifetime could be higher. In Fig.|4l we show relative populations as a function 
of the hold time and derive lifetimes as r ^ 1 s, 0.5 s, and 0.2 s for the n = 3, 4, and 5 MI 
phases, this is shorter than predicted which is possibly due to secondary collisions. For n = 1 
and 2, lifetimes of over 5 s were observed. 

In this work we have used an atomic clock transition to characterize the superfluid - Mott 
insulator transition. This demonstrates the power of applying precision methods of atomic 
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physics to problems of condensed matter physics. Exploiting atomic clock shifts, we could 
distinguish sites according to their occupation number. This allowed us to directly image the 
Mott insulator shells with filling factors n = 1 - 5, and to determined the onsite interaction and 
lifetime as a function of n. 

In the future, this method can be used to measure the number statistics as the system un- 
dergoes the phase transition. One would expect that the spectral peaks for higher occupation 
number become pronounced only at higher lattice depth — an indication of this can be seen 
already in Fig. 1. For low lattice depths, the tunneling rate is still high, but one can suddenly 
increase the lattice depth and freeze in populations fTSl, which can then be probed with high 
resolution spectroscopy. Fluctuations in the atom number could identify the superfluid layers 
between the Mott shells. In addition, by applying a magnetic gradient across the lattice, tomo- 
graphic slices could be selected, combining full 3D resolution with spectral resolution of the 
site occupancy. The addressability of individual shells could be used to create systems with 
only selected occupation numbers (e.g. by removing atoms in other shells). Such a preparation 
could be important for the implementation of quantum gates for which homogenous filling is 
desirable. For atoms other than rubidium, atomic clock shifts are much larger, e.g. for sodium 
by a factor of 30. Therefore, it should be easier to resolve the Mott insulator shells, unless 
the coUisional lifetime of the upper state of the clock transition sets a severe limit to the pulse 
duration lfT9l . 
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Figure 1: 



Fig. l.Two-photon spectroscopy across the superfluid to Mott insulator transition. Spectra for 
3D lattice depths of Erec (open squares), 5 Erec (open diamonds), 10 Erec (open circles), 25 
Erec (solid squares), and 35 Erec (solid circles) are shown. The spectra are offset for clarity. 
The shift in the center of the n=l peak as the lattice depth is increased is due to the differential 
AC Stark shift from the lattice. The dotted lines show gaussian fits of the peaks. 
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Figure 2: 



Fig. 2. Probing the onsite interaction energy, (a) The separation between the n = 1 and n =2 
peaks is shown for lattice depths of V = 25 Erec (squares) and V = 35 Erec (circles). As the 
lattice depth was increased the separation increased from 22(1) Hz to 30(1) Hz. The shaded 
area gives the expected value determined from a band structure calculation, including the un- 
certainty in the scattering lengths. The uncertainty in the measured separation is indicated by 
the size of the points, (b) Location of resonances for all MI phases relative to the n = l phase for 
V = 25 Erec and V = 35 Erec- For low site occupation (n = 1, 2, 3), the separation between the 
resonances is approximately constant, implying constant on-site interaction energy U. For V = 
35 Erec, the separation between the n = 4 and 5 peaks was 22(2) Hz, a 27% decrease from the 
30(1) Hz separation between the n = I and 2 peaks. The slope of the lines is fit to the separation 
between the n = 1 and 2 peaks. 

11 



(a) 



so. 



P4 
II 



<10. 



20. 



(5) 



(4) 



(3) 



I 



(2) 

m 

\ 



(1) 



(b) 



7940 7980 8020 

RF Frequency 670000 Hz] 



6060 





(i) 


(2) _ 


(ii) 


(3) 


(iii) 


(4) 




(5) 



Figure 3: 



Fig. 3. Imaging the shell structure of the Mott insulator, (a) Spectrum of the Mott insulator at 
V = 35 Erec- (b) Absorption images for for decreasing rf frequencies. Images 1, 2, 3, 4 and 5 
were taken on resonance with the peaks shown in (a) and display the spatial distribution of the 
n =1-5 shells. The solid lines shows the predicted contours of the shells. Absorption images 
taken for rf frequencies between the peaks (images i, ii and iv) show a much smaller signal. 
The field of view was 185 fim x 80 fim. 
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Figure 4: 



Fig. 4. Lifetime of individual MI shells. The lifetime for each MI phases can be measured 
independently by adding a hold time before applying the two photon pulse. Spectra are shown 
for hold times of ms (filled circles), 100 ms (filled squares), 400 ms (open circles), 2000 ms 
(open squares). The lattice depth was V = 35 Erec except for the 100 ms hold time, for which it 
was V=34 Erec- The lines show gaussian fits to the peaks, and the spectra were offset for clarity. 
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